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Abstract: We assess the current state of searches at the LHC for additional Higgs bosons 
in light of both direct limits and indirect bounds coming from coupling measurements of the 
Standard Model-like Higgs boson. Given current constraints, we identify and study three LHC 
searches that are critical components of a comprehensive program to investigate extended 
electroweak symmetry breaking sectors: production of a heavy scalar or pseudoscalar with 
decay to tt; bb and tt associated production of a heavy scalar or pseudoscalar with decay 
to invisible final states; and tb associated production of a charged Higgs with decay to tb. 
Systematic experimental searches in these channels would contribute to robust coverage of 
the possible single production modes of additional heavy Higgs bosons. 
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1 Introduction 

Following the discovery of a Standard Model-like Higgs boson at the LHC, the systematic 
search for additional weakly-coupled scalars near the electroweak scale is of paramount im¬ 
portance. A variety of experimental searches have been performed for such extended Higgs 
sectors to date, predominantly targeting new scalars with substantial couplings to the elec¬ 
troweak gauge bosons [1-7] or to down-type fermions [8, 9]. However, these searches are only 
sensitive to a fraction of the interesting parameter space in general extended Higgs sectors. 
This raises the surprising possibility that a large number of additional Higgs bosons have 
been produced at the LHC without leaving signals in existing search channels. 

Moving forward into Run 2 at the LHC, a natural question is how the search for additional 
Higgses should be organized in order to ensure systematic coverage of extended electroweak 
symmetry breaking sectors. Given the proliferation of potential signals, it is useful to con¬ 
sider signatures broadly. Interesting topologies and searches for new Higgs bosons can be 
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classified in terms of simplified models, much in the spirit of simplihed model searches for 
supersymmetry [10-12]. These simplified models may then be combined to provide coverage 
of the parameter space of a given extended Higgs sector. 

A useful hrst step in organizing the search for additional states is to begin with the signa¬ 
tures of a single additional Higgs boson, so that the available decay modes involve Standard 
Model (SM) bosons and fermions, the 125 GeV Higgs, and potentially additional invisible 
decays. In a general extended Higgs sector there may be numerous Higgs bosons beyond the 
SM Higgs. By focusing on one new state at a time, we can characterize the dominant signals 
of an extended Higgs sector if the additional Higgs bosons are well separated in mass, or if the 
additional Higgs bosons are approximately degenerate so that decays between heavy Higgs 
bosons are kinematically disfavored. Having comprehensively covered these signatures, it is 
then possible to systematically expand the picture to consider production and decay processes 
involving more than one heavy Higgs boson. 

Within the space of signatures of a single new Higgs state, further powerful guidance is 
provided by the coupling measurements of the recently-discovered SM-like Higgs boson, which 
constrain the parameter space of extended Higgs sectors. These coupling measurements are 
currently consistent with SM predictions to within the 20—30% level. Such agreement suggests 
that any extension of the Higgs sector must be near an alignment limit in its parameter space, 
wherein the SM-like Higgs boson is closely aligned with the vacuum expectation value (vev) 
that breaks electroweak symmetry and correspondingly exhibits the properties of a SM Higgs 
boson [13-16]. In a given extended Higgs sector, this alignment limit may be approached 
either due to decoupling of additional Higgs states [16, 17], or simply due to the organization 
of dimensionless couplings in the Higgs potential [13-15, 18-20]. Proximity to the alignment 
limit then governs also the couplings of additional Higgs bosons, and may be used as a guide 
to searches for additional Higgses. 

The precise properties of the SM-like Higgs boson and additional Higgs scalars near 
the alignment limit depend on the nature of the extended Higgs sector. The SM Higgs 
boson is a vacuum state in that it carries the quantum numbers of the vacuum. Additional 
neutral Higgs bosons may include pure vacuum states with the quantum numbers of the 
vacuum, allowing mixing with the SM Higgs, but without intrinsic coupling to gauge bosons 
or fermions. Mixing with new pure vacuum states modify the SM-like Higgs boson couplings 
in a model-independent way by diluting all couplings uniformly. The real singlet extension of 
the Higgs sector provides a natural example of a pure vacuum state, with one new CP-even 
Higgs boson whose couplings to SM fermions and gauge bosons are uniformly suppressed in 
the alignment limit. In this case, proximity to the alignment limit implies suppression of all 
production modes for the additional boson at the LHC. 

Alternately, additional neutral Higgs bosons may simply include new vacuum states al¬ 
lowing both mixing with the SM-like Higgs and independent intrinsic couplings to massive 
gauge bosons and fermions. Such vacuum states modify the SM-like Higgs boson couplings 
in a model-dependent way. Such vacuum states arise, for example, in CP-conserving Two- 
Higgs-Doublet Models (2HDM). The physical spectrum of 2HDM includes four additional 
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Higgs bosons - a CP-even scalar H, CP-odd pseudoscalar A, and charged Higgs bosons 
of which the CP-even scalar H without any quantum numbers can be identified as a vac¬ 
uum state. The couplings of these additional Higgs bosons to SM bosons are suppressed in 
the alignment limit, while their couplings to SM fermions are generically unsuppressed but 
depend in detail on the coupling structure of the 2HDM. In this case, proximity to the align¬ 
ment limit implies suppression of production and decay modes involving SM bosons, while 
production and decay via SM fermions may be appreciable. Consequently, typical searches 
such as H ^ ZZ become ineffective. Similarly, H/A —)• rr and —?■ tu searches may be 

effective in some scenarios, but are ineffective whenever the down-type fermionic couplings 
are not substantially enhanced over the SM Yukawas. 

In this paper we articulate a systematic strategy for searching individually for additional 
Higgs scalars in light of the properties of the 125 GeV Higgs. We focus on the phenomenology 
of CP-conserving scenarios with two Higgs doublets satisfying the Glashow-Weinberg (GW) 
condition [21], as this describes the physics of many well-motivated extensions of the Higgs 
sector while still covering many of the key features of models with additional singlets or 
higher electroweak representations. We first summarize the state of limits on 2HDM at the 
LHG from direct searches for additional Higgs states and indirect constraints from Higgs 
coupling measurements (for recent related work, see [14, 22-37]). We then identify and study 
three primary channels where, without being meaningfully constrained by existing searches, 
a second Higgs boson could exhibit 0(1) signals: 

1. The single production of a heavy scalar or pseudoscalar Higgs boson with decay to tt. 

2. The single production of a heavy scalar or pseudoscalar with decay to invisible final 
states. 

3. The tb associated production of a charged Higgs with decay to tb. ^ 

In each case, “single production” includes both resonant production of a single heavy 
Higgs boson and potential associated production modes involving SM fermions in conjunc¬ 
tion with a heavy Higgs boson. The combination of Higgs coupling measurements, ongoing 
searches for heavy Higgses, and the three search channels studied in this work should con¬ 
tribute a rather comprehensive coverage of individual scalar states in extended electroweak 
symmetry breaking sectors. 

The paper is organized as follows: In section 2 we first review the constraints from current 
Higgs coupling measurements on the parameter space of motivated 2HDM, then present the 
combined impact of existing direct searches for heavy Higgs states on the same parameter 
space. In section 3 we consider one of the most pressing signatures of additional Higgs scalars 
in light of current direct and indirect limits: the strong production of a heavy neutral Higgs 

Although the H'^tb process with dileptonic top decay is the subject of a recent CMS search at 8 TeV [38], 
there is room for improvement in the reach, particularly by inclusion of the semileptonic top decay channel 
studied here at 14 TeV. 
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boson followed by decay into tt. This process has a distinctive interference between signal 
and SM tt background [39], but is hampered by the sizable SM background, and further 
complicated by the modest reconstruction resolution available at the LHC of the invariant 
mass of the tt system. Given the challenges of the search for resonant tt production, we also 
briefly consider ti and bb associated production of a heavy neutral Higgs followed by decay to 
ti. In section 4 we turn to invisible decays of heavy Higgses. We consider ti and bb associated 
production of an invisibly-decaying heavy Higgs boson, as proximity to the alignment limit 
renders ineffective traditional searches involving vector bosons. Finally, in section 5 we study 
the tb associated production and decay of charged Higgs bosons. We conclude in section 6 
and reserve details of our Higgs fit, kinematics of the H/A —)• ti process, and top quark 
reconstruction algorithm for a series of appendices. 

2 Direct and Indirect 2HDM Limits 

Direct searches at the LHC and indirect limits arising from Higgs coupling measurements 
impose constraints on the parameter space of 2HDM. As we will discuss in this section, the 
search for additional Higgses is guided by the complementarity of these direct and indirect 
constraints. 

In light of stringent flavor constraints, we focus on (CP-conserving) 2HDM satisfying 
the Glashow-Weinberg condition that all fermions of a given representation receive their 
masses via renormalizable couplings to a single Higgs doublet. There are four distinct possible 
configurations satisfying the GW condition; in this paper we will further focus on the two 
most common, known as Type 1 and Type 2 2HDM. In Type 1 2HDM all SM fermions couple 
to one doublet, while in Type 2 2HDM the up-type quarks and down-type quarks/leptons 
couple to separate doublets. These two types arise most frequently in motivated extensions of 
the SM, including composite Higgs models, little Higgs models, and supersymmetric models; 
the Higgs sector of the MSSM is an instance of the Type 2 2HDM. 

In theories with two Higgs doublets and the most general renormalizable CP- 

conserving potential, there are nine free parameters that remain after minimizing the potential 
and fixing the symmetry breaking vev v‘^ = + = (246 GeV)^. There are various possible 

parameterizations. Here we use the conventions of [14], taking for the free parameters the 
ratio tan j3 = | ($ 2 )1 1 the mixing angle a that diagonalizes the neutral scalar mass matrix, 
the four physical masses {mh,mH,mA,m}{±}, and the dimensionless couplings As^e,?- 

The coupling of the physical states h, H, A, to SM fermions and gauge bosons are fully 
determined by the angles a and /?, while the renormalizable couplings involving three or four 
physical Higgs bosons depend on the additional parameters of the potential. The couplings 
of physical scalars to SM fermions and gauge bosons as a function of a and /3 in Type 1 and 
Type 2 2HDM are summarized in table 1. In this work we will assume that the observed 125 
GeV Higgs is the CP-even scalar h with SM-like Higgs couplings, with the additional Higgs 
scalars iL, A, parametrically heavier. The case of additional scalars lighter than the 125 
GeV Higgs is also quite interesting but qualitatively distinct. 
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Table 1: The coupling of Higgs bosons h,H,A to SM bosons and fermions as a function 
of the angles a and (3, expressed in terms of the alignment parameter Cjs-a = cos(/? — a), 
and = tan (3. The coupling dependence of the charged scalars is the same as the 
pseudo-scalar A. 


It is apparent from table 1 that couplings of the CP-even scalar h become exactly SM-like 
in the limit cos(/3 — a) —>• 0, which coincides with the alignment limit for 2HDM satisfying 
the Glashow-Weinberg condition. In the alignment limit the heavy CP-even neutral Higgs H 
decouples from SM vector bosons, and its couplings become akin to those of the pseudoscalar 
Higgs A. Crucially, the Higgs bosons H, A, and retain couplings to SM fermions in 
the alignment limit. These couplings ensure that the additional states have non-vanishing 
production channels and visible decay signatures involving SM fermions even in the limit 
where the 125 GeV Higgs is exactly SM-like. 

At present, the SM-like nature of the 125 GeV Higgs boson implies proximity to the 
alignment limit commensurate with the precision of Higgs coupling measurements. In order 
to quantify the impact on the (a, (3) parameter space of 2HDM, we perform a global fit to 
recent Higgs measurements reported by the ATLAS and CMS collaborations.^ We provide 
details of our ht procedure in appendix A. In figures 1 we show the result of global hts for 
Type 1 and Type 2 2HDM as a function of tan (3 and cos(/3 — a). We refer the reader to [14] 
for discussion of the physics underlying the shape of these fits. 

The proximity to the alignment limit implied by coupling measurements of the SM-like 

^For this fit and for the interpretation of direct searches for heavy Higgs bosons, we use the programs 
HIGLU/HDECAY [40] to determine the NLO dependence of the h/HjA gluon fusion production cross section 
and partial widths h/H —>■ gg, tt, bh, ss, cc, ng, tt, WW, ZZ, and A ^ gg on the parameters a and /3. We use 
analytic NLO QCD expressions for the partial widths hjH/A —>■ 77 [41] and A —>■ tt, bb, cc [42, 43]. We use the 
program SusHi [44] to determine the NLO bbhjHjA production cross section and validate the HIGLU result for 
gluon fusion. We use MadGraph 5 [45] to determine the LO tth/H/A production cross section with a fe-factor 
of 1.18 [46]. We use leading order results for the partial widths H ^ hh and A — > Zh,TT,gg [47]. 
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Figure 1: Coupling fits in the 2HDM parameter space of tan/3 and cos(/3 — a) in Type 1 
(left) and Type 2 (right) 2HDM. Details of the fit procedure are discussed in appendix A. 


Higgs provides a natural organizing principle for the signatures of additional Higgs bosons. 
The implications for production modes are particularly transparent. In Type 1 2HDM, current 
fits require cos(/3 —a) < 0.4, suggesting that vector associated production modes of H such as 
ZH associated production or vector boson fusion (VBF) are suppressed by at least a factor 
~ 0.2 relative to a SM Higgs of the same mass. In contrast, strong production modes may 
remain appreciable. Gluon fusion production of or A proceeds through fermion loops as 
in the SM, uniformly proportional to cot^/3 in the alignment limit. The same is true of 
tiH/A and bbH/A associated production and tbH^ associated production, indicating that 
these channels remain promising in the alignment limit of Type 1 2HDM. 

In Type 2 2HDM the suppression implied by Higgs coupling fits is even more extreme, 
such that vector associated production modes of H are at most ~ 1% of a SM Higgs of the 
same mass. As in the case of Type 1 2HDM, strong production modes are still appreciable. 
Gluon fusion production of H and A again proceeds through fermion loops, with the top loop 
contribution proportional to cot^ /3 and the bottom loop contribution proportional to tan^ /3 
at leading order in the alignment limit. The tiH/A associated production mode again scales 
as cot^/3, while the bbH/A associated production mode scales as tan^/3. Production of the 
charged Higgs is a function of both tan/3 and cot /3 in the alignment limit. 

The impact on branching ratios of heavy Higgs bosons is somewhat more subtle. As 
discussed in detail in [14], although proximity to the alignment limit implies suppression of 
couplings to SM bosons, these longitudinally-enhanced partial widths are competing only 
with relatively small fermionic partial widths. As such, decays into SM bosons may remain 
appreciable close to the alignment limit. In the exact alignment limit, tree-level decays into 
massive SM bosons (including the 125 GeV Higgs h) vanish in favor of decays into SM fermions 
and the massless gauge bosons.^ 

^We do not consider loop-level decays into massive vector bosons, which are nonzero in the exact alignment 
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Table 2: Hierarchy of heavy Higgs leading LHC production channels that do not vanish in 
the 2HDM alignment limit. 


In table 2 we summarize the the leading LHC production channels for heavy Higgs bosons 
in 2HDM that are non-vanishing in the alignment limit, ordered by their relative size at LHC 
energies. These include resonant production of heavy neutral Higgses by gluon fusion; single 
production of heavy neutral or charged Higgses in association with top and bottom quarks; 
heavy Higgs pair production via Drell-Yan processes; heavy-light Higgs boson production via 
gluon fusion; and heavy Higgs pair production via gluon fusion. Other production modes that 
vanish in the alignment limit are significantly suppressed near the alignment limit, rendering 
them unpromising in the parameter space currently allowed by Higgs coupling hts. We 
likewise summarize the Standard Model decay channels of heavy Higgs bosons in table 3. In 
contrast with production modes, decay modes that vanish near the alignment limit may still 
be appreciable near the alignment limit, given the relatively small partial widths of competing 
decays. 

Given proximity to the alignment limit, there is a natural ordering of searches for ad¬ 
ditional Higgs bosons obtained by combining the dominant production and decay modes. 
Many of the single heavy Higgs boson production channels are covered by existing searches, 
including searches for gluon fusion production of H/A with decay to bb,TT,'y'^, faia as well 
as WW, ZZ, Zh,hh] searches for hbH/A associated production with decay to bb,TT,figL; and 

limit but sufficiently small to avoid influencing the tree-level result. 
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Table 3: Standard Model decay channels of 2HDM heavy Higgs bosons. A checkmark 
indicates that the partial decay width approaches a constant in the alignment limit, while a 
dash indicates that the decay width vanishes in the alignment limit. 


tbH^ associated production with decay to tv and tb. However, several key channels remain 
uncovered, particularly gluon fusion with decay to tt; associated production of bbH/A fol¬ 
lowed by decay to 77 and WW, ZZ, Zh, hh as well as tt] and associated production of tiH/A 
with decay to bb,TT,^^, ^j, as well as WW, ZZ, Zh, hh and tt. Once decay into tt becomes 
kinematically accessible, it becomes one of the primary decay modes of heavy neutral Higgs 
bosons near the alignment limit, and this decay channel may entirely dominate the visible 
signatures of additional Higgses. Similarly, tbH^ associated production with decay to ib is 
likely to be a dominant signature of charged Higgses at the LHC when this decay channel 
is open. Although there is a search for this mode at ^/s = 8 TeV [38], there is room for 
improvement in this channel. 

In addition to decays into SM final states, it is possible for new Higgs bosons to decay 
into non-SM final states. These processes include both invisible decays and potentially visible 
decays that do not fall into the acceptance of existing searches. Given the suppression of vector 
associated production modes in the alignment limit, the most promising potential channels 
are tiH/A or bbH/A associated production with decay to invisible final states. 

To fully characterize the state of coverage by direct searches, we interpret searches 
by the ATLAS and CMS collaborations for heavy Higgs states in the parameter space of 
Type 1 and Type 2 2HDM. The relevant search channels are summarized in table 4. These 
searches present limits in terms of single-channel cross sections times branching ratios that are 
amenable to reinterpretation. Powerful limits on gg ^ H ^ hh and gg ^ A ^ Zh for mod¬ 
erately heavy H,A have also been obtained using multi-lepton and di-photon final states [7], 
but these bounds combine many exclusive channels with non-uniform scaling and acceptance 
across the 2HDM parameter space and cannot be easily reinterpreted in our framework. 
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Channel 

Collaboration 

Reference 

55 ^ ^ 77 

ATLAS, 20.3 fb"! 

[48] 

55 ^ 77 

CMS, 19.7 fb-i 

[49] 

gg ^ ^ ^ TT 

ATLAS, 20.3 fb"! 

[ 8 ] 

66 — ^ TT 

ATLAS, 20.3 fb-^ 

[ 8 ] 

— >■$—>■ TT 

CMS, 19.7 fb-i 

[9] 

66 — > $ — > TT 

CMS, 19.7 fb-i 

[9] 

gg ^ A ^ Zh -A' + (66, tt) 

ATLAS, 20.3 fb"! 

[ 1 ] 

gg —>■ A —>■ Zh —^i£ + bb 

CMS, 19.7 fb-i 

[5] 

gg ^ H ^ hh ^ bb + 

ATLAS, 20 fb-i 

[50] 

QQ —y H —y hh —y hh -\- bb 

CMS, 17.9 fb-i 

[51] 

gg ^ H ^ hh ^ bb + 77 

CMS, 19.7 fb-i 

[52] 

gg ^ H ^ ZZ ^ U 

ATLAS, 20.7 fb-^ 

[53] 

gg —>■ H —>■ ZZ 

CMS, 19.7 fb-i 

[ 6 ] 

gg^ H ^ WW 

CMS, 19.7 fb-i 

[ 6 ] 


Table 4: Relevant ATLAS and CMS searches for heavy Higgs bosons at the ^/s = 8 TeV 
LHC. Here = LT, A. 


For each search, we consider the contribution of Lf or A separately (in contrast to e.g. the 
MSSM interpretation of searches in the rr final state, which includes the sum of contributions 
from h,H, and A). To determine the theory prediction for relevant cross sections times 
branching ratios across the 2HDM parameter space, we obtain the relevant cross sections and 
partial widths as a function of a and j3 as discussed above. Here we assume that the total 
widths of H and A are determined purely by their decays into SM final states. 

In figure 2 we present the state of current direct searches in the exact alignment limit 
cos(/3 — a) = 0 for heavy CP-even neutral scalar H and CP-odd neutral pseudoscalar A as 
a function of tan/3 and in Type 1 and Type 2 2HDM. In the exact alignment limit, 

only production and decay modes involving Higgs couplings to fermions (including gluon 
fusion production and decay into photons arising from top/bottom quark loops) contribute. 
In Type 1 2HDM all production modes involving fermions are suppressed at large tan/3, 
so that existing searches are only effective at low tan /3. The most sensitive search channels 
include inclusive production of H/A followed by decay to 77 , rr. These channels are modestly 
effective near tan/3 = 1 for rriH/A ^ 350 GeV, but lose sensitivity for > 2mt once 

decays into tt go on-shell. In Type 2 2HDM both the gluon fusion and bbH/A associated 
production modes grow at large tan 13, providing additional sensitivity relative to the Type 1 
scenario. Note that the exclusion due to our interpretation of searches in the rr final state is 
somewhat weaker than the comparable MSSM exclusion plot. This is due to the fact that the 
MSSM interpretation combines contributions from h, H, and A, whereas we consider only the 
contribution due to Lf or A individually. In both 2HDM types, the profound weakening of 
limits at low tan /3 in the alignment limit when the H/A —)• tt channel becomes kinematically 
accessible highlights the need for effective searches in the tt final state. 
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H in the Alignment Limit, Type 1 2HDM 



niH (GeV) 


□ ATL □ CMS gg-^O-^yy' □ ATL gg-»cp-»rr 

□ CMS gg-»0^rr □ ATL bb-»0-»rr □ CMS bb^0-»rT 



200 300 400 500 600 


niH (GeV) 

□ ATLgg^C>-»y/y' □ CMS gg-»cb^y/y/ □ ATL gg-»0^rr 

□ CMS gg-»0^rr □ ATL bb-»cb^rr □ CMS bb^0-»rr 

Figure 2: Top: Direct search limits on a heavy 
neutral pseudo-scalar A (right) as a function of 
in a Type 1 2HDM. Bottom: Same as above for 


A in the Alignment Limit, Type 1 2HDM 



ntA (GeV) 

□ ATL gg^<t)-»y/'y' □ CMS gg-»0-»y'y' □ ATL gg-»0-»rr 

□ CMS gg-»<b^Tr □ ATL bb-»<b^Tr □ CMS bb^0->rr 


A in the Alignment Limit, Type 2 2HDM 



rriA (GeV) 

□ ATL gg^O-^y^y' □ CMS gg-»(b^y/'y □ ATL gg-»(b^rr 

□ CMS gg-»0^Tr □ ATL bb-»O^Tr □ CMS bb^®-»rr 


CP even neutral scalar H (left) and CP odd 
mass in the alignment limit cos(/3 — a) = 0 
a Type 2 2HDM. 


As we move away from the exact alignment limit, vector boson associated production 
modes remain unimportant, but decays into vectors can become appreciable. Given the 
sensitivity of searches for heavy scalars decaying into SM bosons, searches in these final 
states become significant relatively close to the alignment limit. In figure 3 we present the 
state of direct searches for H/A with = 300 GeV as a function of tan,0 and cos(/3 — a) 

in Type 1 and Type 2 2HDM. As in the case of the exact alignment limit, for Type 1 2HDM 
sensitivity falls off with increasing tan due to the falling production cross section. The 
strongest limits on H are provided by searches for gluon fusion production of H followed 
by decays into ZZ —)• Ai, although these limits fall off near the alignment limit, where they 
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H Bounds, mH=300 GeV, Type 1 2HDM 



cos(jS-a) 

□ ATL □ CMS gg-»0-»K^' □ ATL gg^<I>-»rT □ CMS gg-»0-»rr 

□ ATL bb-»cp-»rr □ CMS bb-»ct)-»rr □ ATLgg-»H-»hh □ CMS gg-»H-»hh->4b 

□ CMS gg-»H-»hh-»bb// □ ATL gg-»H-»ZZ □ CMS gg-»H^W 



- 0.2 - 0.1 0.0 0.1 0.2 


cos(j8-a) 

□ ATL gg-»<t>-»^'i/ □ CMS gg-^O-^y'y' □ ATLgg^<t>^rr □ CMS gg-»0-»rr 

□ ATL bb-»0-»rr □ CMS bb-»<I>-»rr □ ATL gg-»H->hh □ CMS gg-»H-»hh-»4b 

□ CMSgg^H^hh-»bb// □ ATL gg^H^ZZ □ CMS gg^H^VV 


A Bounds, mA=300 GeV, Type 1 2HDM 



cos(jS-a) 

□ ATL gg-^O-^// □ CMS gg^tlJ-^K/ □ ATL gg-»<t>^rr 

□ CMS gg-»(j5-»Tr □ ATL bb^<I)->rr □ CMS bb-»cp-»Tr 

□ ATL gg-»A-»Zh^llbb □ CMS gg-»A-»Zh^llbb 


A Bounds, mA=300 GeV, Type 2 2HDM 



cos(jS-a) 

□ ATL gg^O-^y;/ □ CMS gg^tl*-)// □ ATL gg-><J>-^rr 

□ CMS gg-»0-»Tr □ ATL bb-»<t>^rr □ CMS bb-»0-»Tr 

□ ATL gg-»A-»Zh^llbb □ CMS gg-»A-»Zh^llbb 


Figure 3: Top: Direct search limits on a 300 GeV CP even neutral scalar H (left) and CP 
odd neutral pseudo-scalar A (right) as a function of cos(/? — a) and tan/3 in a Type 1 2HDM. 
Bottom: Same as above for a Type 2 2HDM. Here we have taken As, 6,7 = 0 in all plots. Note 
the different range of cos(/3 — a) for Type 1 and Type 2 2HDM, motivated by the parameter 
space allowed by Higgs coupling fits. 


are supplanted by searches for H —)• 77 , rr. The strongest limits on A come from searches 
for gluon fusion production of A followed by decay into Zh ^ iibb. These limits likewise 
fall off near the alignment limit, where A —)• 77 , rr provides complementary sensitivity. The 
situation for Type 2 2HDM is comparable to the Type 1 2HDM, save that searches in the 
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H Bounds, mH=500 GeV, Type 1 2HDM 
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□ ATL □ CMS gg-»0-»K^' □ ATL gg^<I>-»rT □ CMS gg-»0-»rr 

□ ATL bb-»cp-»rr □ CMS bb-»ct)-»rr □ ATLgg-»H-»hh □ CMS gg-»H-»hh->4b 

□ CMS gg-»H-»hh-»bb// □ ATL gg-»H-»ZZ □ CMS gg-»H^W 


A Bounds, mA=500 GeV, Type 1 2HDM 



cos(jS-a) 

□ ATL gg-^O-^// □ CMS gg^tlJ-^K/ □ ATL gg-»<t>^rr 

□ CMS gg-»(j5-»Tr □ ATL bb^<I)->rr □ CMS bb-»cp-»Tr 

□ ATL gg-»A-»Zh^llbb □ CMS gg-»A-»Zh^llbb 


H Bounds, mH=500 GeV, Type 2 2HDM 



cos(j8-a) 

□ ATL gg-»<t>-»^'i/ □ CMS gg-^O-^y'y' □ ATLgg^<t>^rr □ CMS gg-»0-»rr 

□ ATL bb-»0-»rr □ CMS bb-»<I>-»rr □ ATL gg-»H->hh □ CMS gg-»H-»hh-»4b 

□ CMSgg^H^hh-»bb// □ ATL gg^H^ZZ □ CMS gg^H^VV 


A Bounds, mA=500 GeV, Type 2 2HDM 



cos(jS-a) 

□ ATL gg^O-^y;/ □ CMS gg^tl*-)// □ ATL gg-><J>-^rr 

□ CMS gg-»0-»Tr □ ATL bb-»<t>^rr □ CMS bb-»0-»Tr 

□ ATL gg-»A-»Zh^llbb □ CMS gg-»A-»Zh^llbb 


Figure 4: Top: Direct search limits on a 500 GeV CP even neutral scalar H (left) and CP 
odd neutral pseudo-scalar A (right) as a function of cos(/? — a) and tan/3 in a Type 1 2HDM. 
Bottom: Same as above for a Type 2 2HDM. Here we have taken As,6,7 = 0 in all plots. Note 
the different range of cos(/3 — a) for Type 1 and Type 2 2HDM, motivated by the parameter 
space allowed by Higgs coupling fits. 


rr final state (either in gluon fusion or bbH/A associated production) become appreciable at 
large tan/3. 

We repeat the process for heavier H/A in Type 1 and Type 2 2HDM with rrijjjj^ = 500 
GeV as a function of tan /3 and cos(/3—a) in figure 4. The limits are generally weaker compared 


- 12 ~ 

















































to mill A = 300 GeV, due both to falling signal cross sections and additional contributions to 
the total width from H/A —)■ tt. However, notable exceptions are the bounds on H ^ hh from 
the CMS —)• /i/i —)• 46 search and the bounds on —)• Zh from the ATLAS A —)• Zh —)• Ubb 
search. In particular, the considerable improvement in H ^ hh ^ Ab sensitivity at high mass 
is due to the boosted kinematics of the 46 final state [51]. 

The combination of diverse searches for heavy Higgs bosons demonstrates considerable 
and complementary coverage across a wide range of 2HDM parameter space, but also high¬ 
lights the substantial holes in existing coverage. In particular, in Type 1 2HDM, searches 
lose effectiveness at large tan /3 due to falling signal cross sections, and more generally lose 
sensitivity near the alignment limit when H/A can decay into ti pairs. In Type 2 2HDM 
there is additional sensitivity at large tan (3 due to enhanced gluon fusion and 66 associated 
production, but sensitivity is poor at moderate tan (3. This is due to a combination of low 
production cross sections and, where kinematically available, missed decays into the tt final 
state. Among other things, these holes demonstrate the need for an effective H/A —?• tt 
search. 

3 Searching for a Neutral Higgs in tt 

As we have discussed, the natural place to look for new Higgs states heavier than about 
350 GeV and with SM-like coupling strength to the top quark is in gg ^ H/A ^ tt. It 
has been known since the seminal work of [39] (and recently emphasized in [37]) that this 
channel provides an interesting and challenging opportunity for hadron colliders. In contrast 
to searches for spin-I or spin-2 ti resonances, the spin-0 signal amplitude interferes with the 
QGD background, producing a characteristic peak-dip structure. As such, existing searches 
for ti resonances cannot be meaningfully reinterpreted to place a constraint on additional 
Higgs bosons in the ti final state. 

In this section we begin by revisiting the analysis of [39] with an eye towards the impact 
of detector effects in a realistic collider environment. Given the size of the SM ti background, 
we introduce a novel technique to efficiently model both detector and event reconstruction 
effects with adequate statistics. We find that because smearing of the reconstructed invariant 
mass of the ti system is typically the same order as (or larger than) the widths and range 
of interference effects of the new Higgs states, the peak-dip structure is largely washed out. 
While the statistical significance of the residual excess can become large at high luminosity, 
systematic effects are likely to render this channel unviable using standard reconstruction 
techniques. 

We are thus motivated to consider ancillary probes of the ti final state that are not subject 
to the same interference effects, namely the associated production channels tiH/A —)• titi and 
bbH/A —?■ bbti. We do not perform complete 14 TeV studies of these channels here, but we 
argue that - based on current 8 TeV trilepton limits - the four top quark channel in particular 
is likely to have sensitivity to moderate-mass scalars. 
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Figure 5: Cross sections vs tt invariant mass for pp — )• <f> — )• tt, where ^ = H (A) in the left 
(right) panel. The dashed black line shows the QCD background, and the different solid lines 
are associated to different values of the mass. 

3.1 pp —7> H/A —)• tt 

We begin by considering the leading-order interference effects between the pp ^ H/A ^ ti 
signal and the SM continuum tt background. In figure 5 we reproduce the differential rates 
for pp ^ H/A —?• ti, combining the parton-level cross sections computed in [39] with the 
parton distribution functions (PDFs) evaluated in [54]. The coupling strengths are set by 
the SM top Yukawa and mt = 173 GeV (the full gg —?• ti differential cross section including 
all interference effects for general 2HDM couplings is given in appendix B). The characteris¬ 
tic peak-dip interference structure is apparent, particularly for heavier (pseudo)scalars; the 
signal-background interference term dominates the pure signal term for all heavy Higgs boson 
masses. This highlights the challenge facing searches for H/A ^ ti at hadron colliders even 
before finite detector resolution is taken into account. 

Given the size of the SM ti background and delicacy of the signal-background interference, 
it is crucial to incorporate detector effects with adequate Monte Carlo statistics. To efficiently 
simulate detector effects, we derive composite smearing functions for ti events as follows; We 
consider seven different reference values for the top quark pair invariant mass and for 
each we generate 10® QCD ti events in Madgraph [45], requiring \m^i; — m/^\ < 0.5 GeV. We 
then shower with PYTHIA6.4 [55] and process the events through Delphes3 [56, 57]. We then 
reconstruct the semi-leptonic ti system using mass-shell constraints as detailed in appendix C, 
thereby obtaining a response function mapping to an distribution. In figure 6 we plot 
histograms of these distributions. Interpolating numerically in and we obtain a 
kernel P(m®j, m^f) against which we can convolve the PDF-smeared parton-level differential 
cross section. This allows us to model the effects of detector resolution and ti reconstruction 
on the peak-dip structure without being limited by Monte Carlo statistics. We plot the results 
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Figure 6: Distribution of after detector effects and tt reconstruction for different values 
of the produced top quark pair mass 
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Figure 7: Cross sections vs tt invariant mass for pp —)• <1> —)• tt, where ^ = H (^) in the 
left (right) panel. Relative to figure 6, we now include detector and reconstruction effects. In 
figure 8 we plot the difference between the background and signal+background curves. 


in the two panels of figure 7 for the scalar and the pseudoscalar. 

Detector resolution and tt reconstruction completely erode the peak-dip structure in the 
presence of a heavy Higgs, leaving behind only modest shifts in the tt invariant mass dis¬ 
tribution relative to the QCD prediction. In figure 8 we plot the difference between the 
smeared invariant mass spectra predicted by QCD with a heavy Higgs boson and pure QCD. 
The best-m^f-bin statistical significances y/Ax^ at 3000 fb~^ and the corresponding S/B are 
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Figure 8: Difference between the background-only and signal-^background cross section 
curves shown in figure 7 for pp —)• <h —)• tt, where ^ = H (A) in the left (right) panel. 
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Figure 9: Left: Best-m^f-bin statistical significance expected at 3000 fb“^ for the scalar case 
as a function of bin width A. Right: the corresponding S/B. Qualitatively similar results 
hold for the pseudoscalar resonance. 


shown for the scalar resonance in hgure 9 as a function of bin size; qualitatively similar results 
hold for the pseudoscalar. From these hgures, we conclude that although the high-luminosity 
LHC will have sufficient statistical power to observe H/A^ti in principle, systematic un¬ 
certainties (even at the percent level) will almost certainly prevent any significant detection. 
Although we have only considered signal and background and leading order (as full next-to- 
leading-order (NLO) expressions for signal-fbackground do not yet exist), it is unlikely that 
the inclusion of NLO effects will significantly alter these conclusions. 

Of course, there is more information in the ti hnal state than just the invariant mass; 
angular distributions and spin correlations may provide additional handles. In appendix B we 
present a parametrization of the ti differential cross section in terms of a well behaved scat- 
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tering variable that affords some additional discrimination between signal and background. 
A full multivariate analysis employing all this ancillary information would increase the sensi¬ 
tivity incrementally, but we do not expect this would substantially alter our conclusions. 

3.2 pp — )> bbH/A — )• bbtt and pp — )• tiH/A — )• titt 

Given the considerable challenges facing a search in the ti hnal state, it is useful to consider 
associated production modes in conjunction with H/A ^ ti. In the alignment limit, vector 
associated production modes for H such as Higgs strahlung or VBF are strongly suppressed, 
while such modes are entirely nonexistent for A. This suggests focusing on fermionic associ¬ 
ated production modes such as tiH/A or bbH/A. The former is appreciable at low tan/3 in 
both 2HDM types, while the latter is appreciable at moderate to large tan /3 in Type 2 2HDM 
such as the Minimal Supersymmetric Standard Model (MSSM). 

Let us first consider ti associated production of H or A followed by decay to ti. On one 
hand, the resulting 4-top final state provides an abundance of promising signal channels. On 
the other hand, the massive three-body kinematics of ti associated production with > 

2mt lead to at most a C>(fb) rate at 8 TeV. Prospects improve significantly at ^/s = 14 TeV, 
but even here the production cross section is at most in the tens of femtobarn. 

There are a variety of searches at ^/s = 8 TeV that are sensitive to the 4-top final state, 
particularly those involving same-sign dileptons (SSDL) or multileptons with additional b- 
tagged jets. To estimate the reach of the 4 top channel at ^/s = 8 TeV, we reinterpret an 
8 TeV CMS multilepton search [58]. We use SRs 8, 18, and 28 of [58]. These signal regions 
are all characterized by trilepton events with a Z-veto, at least four jets with two 6-tags, and 
Ht > 200 GeV, and are distinguished by I/t hr the range 50-100 GeV, 100-200 GeV, and 
> 200 GeV, respectively. We compute the acceptance times efficiency for tiH/A —)■ titi signals 
in these signal regions using Madgraph/Pythia/Delphes as above. Signal regions 8 and 18 are 
the most sensitive, with SR 28 contributing additional sensitivity for larger values of To 

set limits using the observed event counts in [58] , we treat each bin as an independent Poisson 
variable and combine limits from individual bins using a Bayesian algorithm with a flat prior 
on signal strength, marginalizing over a normally-distributed background uncertainty. We 
neglect potential uncertainties on the signal cross section. 

The 8 TeV data is insufficient to set a limit on a SM-like Hti coupling. For example, we 
find for mn = 350 GeV that the combination of SRs 8, 18, 28 exclude cj • Br > 160 fb. By 
contrast, for mn = 350 GeV and tan/3 = 1 we have a{pp —)■ tiH) ~ 5 fb at ^/s = 8 TeV, 
placing the signal cross section more than an order of magnitude below the current exclusion 
in this channel. While the inclusion of additional channels such as SSDL would improve 
sensitivity, the disparity between signal cross section and exclusion too great to hope for 
meaningful sensitivity in this channel at ^/s = 8 TeV. Although we have explicitly considered 
the case of pp —)• tiH, the rate for pp tiA is comparable, and degenerate H/A in the 
alignment limit would lead to a doubling of the signal. 

At y/s = 14 TeV the prospects for a search in the 4-top channel improve considerably, 
as the tiH/A associated production cross section increases by an order of magnitude for 
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Figure 10: Left: Projection for the excluded cross section times branching ratio for pp —)• 
ttH —7- tttt at y/s = 14 TeV, 3000 fb“^ as a function of rriH- The projection is obtained from 
the naive scaling of the combination of + 2b channels (solid blue line) and with a further 
\/2 improvement from the addition of SSDL channels (dashed blue line), and is compared 
to the signal cross section a{pp —)• ttH) at tan/3 = 1 (solid black line). Right: Projected 
exclusion for a CP-even neutral scalar H in the alignment limit of a Type 1 2HDM from 
the naive scaling of the combination of “ii + 2b channels (solid blue) and with a further \f2 
improvement from the addition of SSDL channels (dashed blue). 


mil/A ^ relative to \/s = 8 TeV. However, reliably estimating sensitivity at the level 
of a theory study is challenging since the largest backgrounds to SSDL and multi-lepton 
searches in this final state typically originate from tt or WjZ -|- jets events with an additional 
fake lepton. As such, we estimate the 14 TeV uxBr exclusion reach with luminosity and 
background cross section rescaling of [58], assuming efficiencies remain comparable to 8 TeV. 
In particular, we rescale background cross sections by the ratio of tt cross sections at 14 TeV 
and 8 TeV, since tt + lepton fakes comprise the dominant background in signal regions 8, 18, 
and 28 of [58]. 

The impact of this projected sensitivity at y/s = 14 TeV is shown in figure 10, both for 
the naive scaling of the sensitivity from the combination of 3i -I- 2b channels, and with an 
additional factor-of-\/2 improvement in the a-Br reach to emulate the potential improvement 
from including SSDL channels. Even at ^/s = 14 TeV this remains a challenging channel, but 
offers hope for meaningful sensitivity at low tan (3 for 2mt < mn < 500 GeV. 

We turn next to bb associated production of Lf or A followed by decay to tt. This process 
may be significant in Type 2 2HDM where the bhH/A associated productiion grows with 
tan/3. However, in this case the partial width H/A^ti also falls as with tan/3, suggesting 
the rate for bbH/A —)• bbtt will peak at moderate values of tan/3. To estimate the LHC 
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Figure 11: Left: Projection for the excluded cross section times branching ratio in femto- 
barns for pp —?■ bbH —)• bbtt at ^/s = 14 TeV, 3000 fb“^ as a function of mn (solid blue line). 
For comparison we show the signal cross section a{pp —)• bbH) at tan/1 = 1 (solid black line). 
Right: Contours of a ■ Br(pp —)• bhH —)■ bbtt) in femtobarns for a CP-even neutral scalar H in 
the alignment limit of a Type 2 2HDM. 


sensitivity to pp —)■ bhH/A —)• bhti, we design and simulate a search in the semi-leptonic final 
state at ^/s = 14 TeV. In contrast to the multi-lepton search for the four-top final state, the 
dominant backgrounds for a ttbb final state can be reliably simulated with available Monte 
Carlo techniques. 

We generate parton level signal and backgrounds events using MadCraphS [45] to leading 
order with CTEQ6L1 PDFs [59]. The events are showered with PYTHIA6.4 [55] and DelphesS 
[56, 57] is used to simulate detector effects. Jets are reconstructed using the anti-ZcT algorithm 
with R = 0.5 and are required to satisfy pUp > 20GeV, jr/-^] < 4.5. Charged leptons (electrons 
and muons) are required to have p^ > 15GeV, jr/^j < 2.5, (Ai? = 0.3) < 0.1. The 

6 -tagging efficiency is chosen to be 70% with a 25% (2%) mistagging rate for charm (light) 
jets [60]. 6-tagged jets satisfy > 40GeV, \7f\ < 2.5. 

In addition to the single lepton requirement, we require at least 6 jets with at least 4 
6 -tags in the hnal state to suppress SM backgrounds. We also apply a missing transverse 
energy cut of I^t > 30 GeV and veto events with more than one charged lepton. Top quarks 
and W bosons are reconstructed from the mass-shell constraints, with small corrections for 
detector effects (for details, see appendix C). After top quark reconstruction, we require 
that the signal events satisfy < 5.0, where x is a variable characterizing the quality of the 
reconstruction (see appendix C). The irreducible background is pp —)• ttbb, while the dominant 
reducible backgrounds for this analysis are pp —)• ttbj and pp —)• tijj, with light jets faking 
bottom quarks. The backgrounds from tth, ttZ, single top production and vector boson plus 


19 





multijets are subdominant [61]. 

To set an exclusion limit, we use the likelihood function 



(3.1) 


where Xj is the binned rriti distribution predicted by the model (with or without signal) and 
Uj is the observed distribution. The 2a exclusion bound is obtained [62] from 



(3.2) 


The results are shown in figure 11. The excluded cross section ranges from 200-70 fb for 
350 GeV < mu ^ 700 GeV, while a{pp —)• bbH) < 10 fb for tan /3 = 1 across the same range. 
Although the production cross section grows with tan/3, the branching ratio to ti falls, so 


that the peak rate a{pp —)• bbH —)• bbtt) ~ 50 fb is obtained around tan/3 ~ 5. Based on 


the results of our preliminary simulation, a meaningful limit cannot be set in Type 2 2HDM. 
However, given that sensitivity is of the same order as the peak rate for 350 GeV < mn ^ 500 
GeV, this channel deserves further experimental study at 14 TeV. 

Finally, we note a third associated production channel that may prove useful in the 
hunt for H/A ^ tt, although we do not study it in detail here. The rate for electroweak 
production of the single-top t{q)H/A final state via a t-channel W boson exceeds that of 
tiH/A production around ra^ijA — 340 GeV [63] and has the same parametric scaling as a 
function of a and /3. Gonsequently, this suggests the rate for pp —)■ t{q)H/A —)■ t{q)tt exceeds 
that of pp —)• tiH/A —?■ titi in the entire region of interest for H/A ^ ti. The resulting 
three-top final state is particularly amenable to a search for same-sign dileptons with two or 
more 6-tagged jets, and may provide a complementary probe of H/A ^ ti at y/s = 14 TeV.^ 

4 Searching for an Invisible Neutral Higgs 

In addition to decays of new scalars to SM fermions, it is also interesting to consider invisible 
decays, which have been actively studied for the SM Higgs boson following Run 1 (with current 
upper limits around 30% coming from the VBF channel [65]). As with the SM Higgs, any 
observation of an invisible width for a second Higgs state would provide a window into new 
physics, possibly signaling the first laboratory production of dark matter. Invisible decays 
could also provide a discovery mode for new Higgses, since SM backgrounds can be strongly 
suppressed with a large missing energy cut. In this section, we continue the study of new 
scalars in channels involving tops and bottoms, adding the ingredient of large missing energy. 

In section 4.1 we briefly discuss some of the theoretical motivation for searching for 
invisibly decaying new scalars in association with tops and bottoms. There are many possible 

'^For a recent study of single top production in association with the Higgs state at 125 GeV, see [64]. 
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UV completions. For simplicity we focus on one amusing example, the MSSM Higgs sector 
benchmark points often used in reporting limits on H/A —)■ rr. Due to choice of neutralino 
mass parameters in these benchmark points, the new neutral scalars associated with the 
second Higgs doublet possess branching fractions into pairs of the lightest i?-odd particle. 
With some variation of the parameters, the invisible branchings can become substantial. 

In section 4.2, we place a limit on bbH —)• 66 + by reinterpreting an 8 TeV sbottom 
search [66]. We give the corresponding limit in the parameter space of Type 2 2HDM, where 
the bbH coupling is tan/3-enhanced, and argue that it is likely to be a stronger limit on this 
parameter space than one derived from monojet searches. Previously, in a study focused on 
dark matter simplified models. Ref. [67] obtained a limit on the 66-associated invisible scalar 
channel by reinterpreting an ATLAS effective operator study at 8 TeV [68] . We have checked 
that the sbottom and effective operator cut flows provide very similar reach, and present our 
results on natural parameter spaces for new scalar searches. 

In section 4.3 we study the reach of semileptonic ttH —)• ft + a-t 14 TeV. Previously, 

Ref. [69] obtained a limit on this channel at 8 TeV by reinterpreting a CMS stop search [70]. 
Ref. [71] also performed an 8 TeV analysis, and furthermore estimated a 14 TeV limit by a 
parton-level reinterpretation of an ATLAS stop study [72]. We complete the phenomenological 
analysis of semileptonic ttH tt + performing a full 14 TeV analysis with optimized 
cuts and detector simulation. We also argue that this channel is likely to be competitive with 
monojet searches for new invisibly-decaying scalars with masses below 2mt. 

4.1 Models with H/A^ 

In the most model-independent spirit of simplified models, any search for an invisibly de¬ 
caying new scalar is of interest: the topologies are simple and capitalize on the small SM 
backgrounds. The production channels studied here are motivated by the alignment limits of 
new weakly-coupled scalar models, which preserve the SM-like properties of the light observed 
Higgs boson at the expense of suppressing traditional invisible scalar searches involving gauge 
couplings such as ZH ii + and qqH qq + R is not difficult to add additional 
theoretical structure, such as dark matter candidates or particles that are long-lived on detec¬ 
tor timescales, into which new scalar states may decay invisibly with substantial branching 
fraction. 

Searches in the ti-\-I^T channel are most effective in models where new scalars couple to ti 
with coupling y[ ~ In the alignment limits of Type 1 and Type 2 2HDM, = cos /3 

and y't/yt^ = cot/3, respectively, so the LHC reach is strongest when tan/3 ~ 0(1). A simple 
model of invisible decays can be obtained by coupling the second doublet $2 to a new massive 
singlet scalar through a portal-type coupling: 

VZiK\^2\'^S^ + ]^yilS^ + ... (4.1) 

In the alignment limit, the dominant decays of the neutral components H, A of <I >2 will be 
into Standard Model fermions and S pairs. Let us assume that the singlet mass is small and 
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that mu ^ 2mt, so that the dominant SM decays are into bb. Then the ratio of invisible to 
visible partial widths is approximately: 

_ kVcos^/3 _ 

^bb ^’rnjjy'b^ 

In Type 1 2HDM, = y®'^cos/3 in the alignment limit, so this ratio is 0(1) already for 
K ~ 0.1 in the range of niH considered. In Type 2 2HDM, is tan/3—enhanced in the 
alignment limit, but for tan/3 ~ 0(1) the invisible decay is still substantial for small k. 
Simple models of this type are most effectively probed by the tt + channel. 

Since yf^ is small, searches for bb + I^t are most effective in cases where the new scalar 
has enhanced coupling to bb. The most well-known example is the tan /3 enhancement of the 
Type 2 2HDM. In the toy model above, the ratio of partial widths is suppressed by (tan/3)“^, 
so K must be large in order to obtain a substantial invisible width. Another Type 2 2HDM 
with the potential for invisible decays is the MSSM. In fact, traditional benchmark scenarios 
used to study H/A^tt actually have regions of parameter space where H/A ^ inv occurs 
at non-negligible rates [73]. 

Supersymmetric Higgs bosons couple to neutralinos in the form Higgs - Higgsino - Elec- 
troweakino. Therefore, for y, ~ Mi <C (i = 1 or 2), the heavy neutral Higgs states can 
decay invisibly into the lightest neutralino through its gaugino/higgsino mixings. The invisi¬ 
ble branching ratios for H are maximized for tan (3 ~ 5; at higher values H ^ bb dominates, 
while at lower values we deviate substantially from the alignment limit and H —?■ VV becomes 
important. In figure 16 we show the H/A ^ inv branching ratios in an MSSM benchmark 
point with tan^d = 5, M 2 = 300 GeV, Mi = (5s^/3c^)M2 = 143 GeV, and decoupled gluino 
and scalars. Sizable invisible branching fractions are possible for both states when the LSP 
is well-mixed and the ti channel is kinematically forbidden. 

We do not pursue model building further in this work, but for completeness we note that 
LUX limits [74] rule out most of the interesting parameter space if the neutralino is dark 

matter and y > 0. If it is only a subcomponent of dark matter, or if it is stable on detector 

timescales but decays outside of the detector (say, through RPV couplings, or to a gravitino), 
then the direct detection limits do not apply. Another intriguing possibility is that y < 0. 
Ref. [75] observed that for \y\ ~ rnisp-, and tan/3 ~few, there is a blind spot on the mu axis 
where the tree-level direct detection amplitudes from SM Higgs and MSSM Higgs exchange 
cancel with each other. 

4.2 bb + l^T at 8 TeV 

In this section, we reinterpret the ATLAS sbottom search (pp bb + I^t) performed in 

Ref. [66] into a limit on pp —?■ Hbb, where H is a new heavy neutral CP-even scalar that 

decays invisibly. We unfold the signal efficiency factors from the cross section limits cTvis 
given in Table 7 of [66]. We focus on Signal Region A (SRA), defined by the cut flow of Table 
1 of [66], because in that region the leading two jets are b-tagged. SRA is divided into hve 
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Figure 12: Contours of the H/A ^ inv branching ratios in an MSSM benchmark point with 
tan/3 = 5, M 2 = 300 GeV, Mi = (5s^/3c^)M2 = 143 GeV, and decoupled gluino and scalars. 
The gray region denotes the LEP exclusion on light charginos. 
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Figure 13: Left: the 8 TeV uxBr limit on the production of invisibly decaying new scalars 
in association with bottom quarks. Right: the corresponding limit on the Type 2 2HDM 
plane. The lighter blue region in the upper left is the bb + exclusion and the darker blue 
region is the estimated monojet exclusion, rescaling the results of [76]. 


subregions based on the contransverse mass cut [77] (with ISR correction given in [78]), and 
the ATLAS observed limits on the cross sections in SRA range from 0.26 — 1.9 fb depending 
on the mcT cut. 

For signal generation and detector simulation we use Madgraph/Pythia/Delphes matched 
to one jet. We optimize the unfolded (ypp^hbH ^ —)• inv) over the mcT cut, finding that 

the softest cut used in the sbottom search {mcT > 150 GeV) gives the best limit on bbH 
until iriH ^ 600 GeV. We attempt to validate our analysis by reproducing the sbottom 
signal efficiencies given in the ancillary data of [66]. We find that near the exclusion limit 
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our analysis yields signal efficiencies approximately a factor of 2 better than those found by 
ATLAS. Therefore, to be conservative we assign a factor of (1/2,2) uncertainty band to our 
bbH limit. 

Our results for the cr x Br limit are given in the left-hand panel of figure 13. The cross 
sections are relatively large, suggesting that this search mode is most effective in constraining 
new scalar models where the scalar coupling to bottom quarks is enhanced over that of the 
SM Higgs. Such enhancements occur, for example, in the alignment limit of Type 2 2HDM, 
where at large tan the Hbb and Abb couplings are a factor of tan /3 larger than the SM Higgs 
hbb coupling. In the right-hand panel of figure 13, we reinterpret the cj x Br limit on the 
parameter space of Type 2 models with Br(iL —)■ inv) = 1. 

We note that a monojet signal arises in this scenario by closing the bottom loop and 
radiating an additional jet. In the case of + the monojet signature obtained in this way 
is competitive and can outperform reinterpreted stop searches [71]. However, in the bb + 
case, we expect that monojet is less powerful for two reasons. First, unlike the ti + case, 
the bottom quarks in the final state do not suffer a large phase-space suppression. Second, 
closing the loop costs a factor of mj in the amplitude, which is a large suppression in the bh 
case, even if the coupling to bottom quarks is tan /3-enhanced. (Moreover, in some models, 
like Type 2 2HDM, the new scalar coupling to top quarks is tan/3-suppressed.) In the dark 
blue region of the right-hand panel of figure 13 we estimate the monojet exclusion by rescaling 
the results of [76] into the Type 2 plane, and we see that it is much weaker than bb + I^t- 
Our rescaling uses the ratio of the LO T{H —)• gg) loop functions in the 2HDM relative to 
the SM, and is therefore rather crude. However, it is likely to be conservative, since the pT 
spectrum of the extra jet in the bottom-dominated process is harder than that of the jet in 
the top-dominated case [79], suggesting that the actual monojet limits will be weaker. 

4.3 ti + at 14 TeV 

We now turn to invisible Higgs states produced in association with top quarks, where limits 
at 8 TeV have been obtained [69, 71] by reinterpretation of a CMS stop search [70], and 
perform an optimized projection for the reach of the 14 TeV LHC. We simulate signal and 
background at leading order in Madgraph/Pythia/Delphes, and apply cuts requiring large 
missing energy, four jets including two 6-tags, and a veto on n/ep 7 ^ 1- The dominant SM 
background processes are semileptonic and dileptonic ti, Zti —)• uuivjjbb, and Wjjbb. 

For our preselection and jet selection cuts, we require I^t > 250 GeV, 1 lepton, at least 
four central jets with pT > {130, 50, 50,30}, respectively, two 6-tags, and A(/)(j, I^t) > 0.8 for 
the two hardest jets. In stop searches, hard cuts on the transverse mass rriT and the variable 
[80] may be used to suppress semileptonic and dileptonic ti, respectively [70]. The 
distributions of I^t, w-t, and after the jet selection cuts are given in Figs. 14, 15, and 16 
for backgrounds and signal for two values of the scalar mass. Subsequently we optimize cuts 
on I^T, mr, and choosing > 300 GeV, rriT > 140 GeV, and > 200 GeV. After 
all cuts, Zti is the dominant background. Subsequently we validate the Zti backgrounds 
against the stop search study in Ref. [70]. 


- 24 


mH=200GeV mH=500GeV 




Figure 14: The I^t distributions of the signal and backgrounds after the jet selection cuts 
described in the text. Left: niH = 200 GeV. Right: ttih = 500 GeV. 




Figure 15: The rriT distributions of the signal and backgrounds after the jet selection cuts 
described in the text. Left: niH = 200 GeV. Right: tuh = 500 GeV. 

Our projected limit {CLg = 0.05) with 300fb“^ is given in the left-hand panel of figure 17. 
It is also interesting to compare with the monojet reach, which we take from the scalar analysis 
performed in [76]. Unlike in the bb+I^T case, the U+I^t search is most effective when the new 
states have SM-like couplings to top quarks. Therefore, rather than comparing the reaches in 
the parameter space of a Type 2 2HDM, we compare the reaches relative to a SM-like fiducial 
model with a new scalar decaying invisibly with unit branching ratio. (Of course, above the 
ti threshold, such a large invisible branching ratio would require a very large coupling to the 
invisible states; the fiducial model is meant only for the comparison of the two search modes.) 
The right-hand panel of figure 17 shows that the associated production mode is expected to be 
competitive with monojet limit until around the tt threshold, where the threshold generates 
a bump in the SM monojet cross section. 
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Figure 16: The distributions of the signal and backgrounds after the jet selection cuts 
described in the text. Left: mn = 200 GeV. Right: mn = 500 GeV. 




Figure 17: Left: the 14 TeV, 300 fb“^ expected limit on the crxBr of invisibly decaying 
new scalars in association with top quarks. Right: comparison with the expected limit from 
monojet exclusion [76], both normalized to the respective SM cross sections. 


5 Searching for a Charged Higgs in tb 

Thus far we have focused largely on the signatures of the vacuum states H and A near the 
alignment limit. In this section, we turn to the signatures of charged Higgs bosons in the 
alignment limit, which may provide an alternative handle on 2HDM in regions of parameter 
space where H and A are hard to find. In particular, we analyze the LHC reach for a charged 
scalar that couples to the SM top and bottom quark through a Yukawa interaction of 
the form 

Aff = VtbH^i {Pl sin 9 + Pr cos 9)9 + h.c. (5.1) 

Near the alignment limit, ih associated production with decay to tb is the dominant channel 
for single production of a charged Higgs boson, as the Wh mode vanishes in the alignment 
limit and decays into tb swamp those into tv. As such, we focus on the process pp —>■ 
H^ib{H~tb) + X with > tb{H~ —)• tb), and we employ the semi-leptonic decay of the tt 
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pair.^ The CMS collaboration recently published a search for the charged Higgs at 8 TeV via 
the same production channel, using the dileptonic decay mode of the top pair [38]; our aim 
is to forecast sensitivity at y/s = 14 TeV and demonstrate the added sensitivity available in 
a semi-leptonic search. Our results are insensitive to the value of 6 in Eq. (5.1), so for the 
results shown here we set 0 = 0®. 

To suppress SM backgrounds, we require at least 4 6-tagged jets in the final state. The 
the irreducible background is 


pp —^ ttbb, (5.2) 

while the dominant reducible backgrounds are 

pp tibj, (5.3) 

pp tijj, (5.4) 

with light jets faking bottom quarks. The tih, ttZ, single top production and vector boson 
with multijets backgrounds are comparatively negligible [61]. 

We generate parton level signal and backgrounds events using MadGraph5 [45] to leading 
order with CTEQ6L1 pdfs [59]. The events are showered with PYTHIA6.4 [55] and Delphes3 
[56, 57] is used to simulate the detector.^ We require at least 6 jets with at least 4 6 -tags, and 
require the leading 6 -jet to have pT > 150 GeV. We also apply a missing transverse energy 
cut of I^T > 30 GeV and veto events with more than one charged lepton. 

Top quarks and W bosons are reconstructed from the mass-shell constraints, with small 
corrections for detector effects (for details, see appendix G). After top quark reconstruction, 
we require that the signal events satisfy < 5.0 and > 0.9, where y is a variable 

characterizing the quality of the reconstruction (see appendix G), and 6 i ( 62 ) is the leading 
(sub-leading) 6 -jet which is not recognized as emerging from a top quark decay. 

The charged scalar invariant mass is reconstructed from the leading 6 -jet and the leading 
reconstructed top quark. In figure 18, we show the tb invariant mass distribution of the 
backgrounds and the signal from a 700 GeV with ytt = 1 at 14 TeV LHG with 3000 fb“^ 
integrated luminosity. 

To obtain the exclusion and discovery bound, we again use the likelihood function given 
by Eq. (3.1), where now where Xj is the binned mtb distribution predicted by the model (with 
or without signal) and rij is the observed distribution. The 2cr exclusion bound is obtained 

®Early investigation of this channel at the LHC can be found in Ref. [81]. 

®The possibility of using this channel to investigate the chirality structure of the H^tb vertex has been 
studied in [82-85]. 

^Jets are reconstructed using the anti-fey algorithm with R — 0.5 and are required to satisfy > 
20GeV, \rf\ < 4.5. Charged leptons (electrons and muons) are required to have py > 15GeV, |p^| < 
2.5, 7iso,y (Ai? = 0.3) < 0.1. The 6-tagging efficiency is chosen to be 70% with a 25% (2%) mistagging 
rate for charm (light) jets[60]. 6-tagged jets satisfy py > 40GeV, |p^j < 2.5. 
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Figure 18: The tb invariant mass distribution of the backgrounds and the signal from a 700 
GeV with yti, = 1 at 14 TeV LHC with 3000 fb“^ integrated luminosity. 



Figure 19: The 2a exclusion and the 5cr discovery bound of the charged Higgs via searching 
the tb resonance in the ttbb channel at 14 TeV LHC with 300 fb~^ (left panel) and 3000 fb~^ 
(right panel) integrated luminosity. 


as in Eq. (3.2), while the 5a discovery reach is obtained from 



L{b\ys + b) \ 
L (/is + b\iJ,s + b)) 


= 5. 


(5.5) 


In figure 19, we show discovery and exclusion curves for the coupling constant ytb as a function 
of niH- We have checked that these results are insensitive to the 9 angle in Eq. (5.1). In the 
Type 2 2HDM, Eq. (5.1) can be written as 


^eS 



V 


t {Pimt cot /3 + PRirib tan /3) 6 + h.c. 


(5.6) 


As such, the constraint on ytb shown in figure 19 is translated into a constraint to tan f3 in 
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figure 20. 



Figure 20: The 2(j exclusion (left panel) and the 5 (t discovery (right panel) bound of the 
charged Higgs via searching the tb resonance in the ttbb channel at 14 TeV LHC with 300 
fb“^ and 3000 fb“^ integrated luminosity. 


While there remains a hole in coverage at moderate values of tan/l, there is nonetheless 
considerable sensitivity for heavy charged Higgses in this channel at ^/s = 14 TeV. Note 
also that the reach of the semi-leptonic search at ^/s = 14 TeV is comparable to the naive 
extrapolation of the ^/s = 8 TeV CMS di-leptonic search, suggesting that an optimized search 
for charged Higgses can effectively employ both semi-leptonic and di-leptonic final states to 
constrain pp —>■ H^tb{H~tb) -|- X with —)■ tb{H~ —)• tb). 


6 Conclusions 

The hunt for the rest of the Higgs bosons is entering a new phase, as an ever-broadening 
set of direct searches at the LHC begins to constrain the parameter space of extended Higgs 
sectors. In this work we have attempted to identify and analyze some of the most promising 
open channels in existing coverage of heavy Higgs bosons consistent with properties of the 
observed SM-like Higgs. These channels are the production of a heavy scalar or pseudoscalar 
with decay to tt] bb and tt and associated production of a heavy scalar or pseudoscalar with 
decay to invisible final states; and tb associated production of a charged Higgs with decay to 
ib. 

Heavy scalars or pseudoscalars decaying into tt constitute a significant gap in existing 
coverage of extended electroweak symmetry breaking scenarios. Taking into account the 
effects of detector resolution and tt reconstruction, we have found that searches for resonant 
production of heavy Higgses with decay into tt are likely to be systematics-limited at the LHC. 
We have correspondingly proposed several ancillary channels involving associated production 
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that may provide complementary sensitivity. The most promising is ttH/A production with 
H/A ti. 

Searches involving missing energy provide an effective probe of intriguing scenarios where 
a heavy Higgs decays invisibly. We demonstrated that bhH/A and tiH/A are valuable pro¬ 
duction channels in which to search for H/A ^ inv. Furthermore, their reach is expected to 
be competitive with ~ or better than - monojet searches in some models and mass ranges. 

Finally, heavy charged Higgses in the alignment limit decay dominantly to tb if this 
channel is open, and the natural strong production channel for charged Higgses with mfj± > 
mt + rrih is in association with tb. We studied the reach in the semileptonic channel, where the 
system can be completely reconstructed, and hnd considerable sensitivity to heavy charged 
Higgses that can complement existing searches in the dileptonic channel. 

In conjunction with precision Higgs coupling measurements and existing direct searches 
for heavy Higgs bosons, these searches can maximize the LHC discovery potential for the 
most well-motivated extensions of the Higgs sector. 
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A Higgs Couplings Fit 


We consider only a single SM-like Higgs boson of mass nih = 125 GeV whose couplings to SM 
fermions and gauge bosons are modihed relative to an SM Higgs boson of the same mass by 
coupling modifiers Ki = gi/gi,sM, where i = t,b, r, W, Z, g, 7 . We treat the loop-induced cou¬ 
plings to gluons and photons independently from variations in the fermion couplings to allow 
for new degrees of freedom running in loops. Eor simplicity we assume custodial symmetry 
so that Kz = K-w- We do not consider a potential invisible width or couplings with non-SM 
tensor structure. 

We construct likelihoods for a Higgs coupling fit using data from Higgs analyses reported 
by the ATLAS and CMS collaborations. Single-channel likelihoods are constructed for each 
Higgs analysis using a two-sided Gaussian where the central value corresponds to the best-fit 
signal strength modifier g reported in the analysis and the variance on each side corresponds 
to the Icj error on the signal strength modifier, i.e. 


Cf{g) oc exp 


' -{g- inf ' 

- 2(uf)2 _ • 


(A.l) 


This two-sided likelihood accommodates the often sizable non-gaussianities found in low- 
statistics channels. 
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The theory prediction for the signal strength modifier is constructed by summing over 
the production and decay modes considered in the analysis (each of which is a function of 
the coupling modifiers Kj), weighted by the relative contribution e of each production mode 
to the analysis. These relative contributions are extracted from experimental publications or 
inferred from the literature where appropriate. We neglect uncertainties on the values of e. 
We consider experimental analyses for which a single decay mode dominates the analysis, so 
that the signal strength modifier for a single experimental channel is given by 


= 



BR 

BRsm' 


(A.2) 


where the index a runs over the gluon fusion, vector boson fusion &: associated vector pro¬ 
duction, and associated tt production modes. The set of ATLAS and CMS Higgs analyses 
used to construct our coupling fit (with corresponding best-fit signal strength modifiers, la 
errors, and relative efficiencies) are enumerated in tables 5 and 6. 


Channel 


Aa 

ea (GGH, VBF/VH, TTH) 

Vh bb (Oi) [86] 

7/8 TeV 

^•«jO_o,52 

(0.0,1.0,0.0) 

Vh bb (H) [86] 

7/8 TeV 

1 1 7 + 0.66 
-‘^•-‘^'- 0.60 

(0.0,1.0,0.0) 

Vh bb (2£) [86] 

7/8 TeV 

0.94t°;«8 

(0.0,1.0,0.0) 

tth — >• bb [87] 

7/8 TeV 


(0.0,0.0,1.0) 

h^TT (jj) [88] 

7/8 TeV 

3.6t?;° 

(0.60,0.4,0.0) 

h^TT (G) [88] 

7/8 TeV 

U - y _ 0.9 

(0.65,0.35,0.0) 

h ^ TT {ii) [88] 

7/8 TeV 

3-0tn7 

(0.65,0.35,0.0) 

hjj TT (jj) [88] 

7/8 TeV 

1 4 + 0-9 
^■+- 0.7 

(0.15,0.85,0.0) 

hjj TT {£j) [88] 

7/8 TeV 

1 0+°® 

(0.12,0.88,0.0) 

hjj TT (U) [88] 

7/8 TeV 

1 S+1-1 

(0.10,0.90,0.0) 

h WW (Oj) [89] 

7/8 TeV 


(0.98,0.02,0.0) 

h WW (Ij) [89] 

7/8 TeV 

0.96tHo 

(0.87,0.13,0.0) 

h WW (2j ggH) [89] 

7/8 TeV 

-L-^U_q 34 

(0.75,0.25,0.0) 

h WW (2j VBF) [89] 

7/8 TeV 


(0.13,0.87,0.0) 

h^ZZ (ggH) [90] 

7/8 TeV 

1.66t°i 

(1.0,0.0,0.0) 

h^ ZZ (VBF-tVH)[90] 

7/8 TeV 

0.26tl:l 

(0.0,1.0,0.0) 

h^-f-f (ggH) [91] 

7/8 TeV 

1 00+0.38 
-i--oz_o 38 

(1.0,0.0,0.0) 

h ^ 77 (VBF) [91] 

7/8 TeV 

0.8t°;^ 

(0.0,1.0,0.0) 

h ^ 77 (WH) [91] 

7/8 TeV 


(0.0,1.0,0.0) 

h ^ 77 (ZH) [91] 

7/8 TeV 

^•-‘-- 0.1 

(0.0,1.0,0.0) 

h — >■ 77 (ttH) [91] 

7/8 TeV 

1 6+2-^ 
■*-'^-1.8 

(0.0,0.0,1.0) 


Table 5: ATLAS Higgs analyses used in constructing coupling fits. The best-fit signal 
strength modifier is denoted by jl with corresponding Tlu errors. The relative contributions 
e are reported for production initiated by gluons via gluon fusion (GGH), weak gauge bosons 
via vector boson fusion or vector associated production (VBF/VH), and top quarks via tt 
associated production (TTH). 


- 31 - 










Channel 


Aa 

ea(GGH, VBF/VH, TTH) 

h ^ bb [92] 

7/8 TeV 

1 n+0.53 
-*-■'^-0.50 

(0.0,1.0,0.0) 

tth —>■ bb [93] 

7/8 TeV 

0.67;1;35 

(0.0,0.0,1.0) 

h-)-TT (0,lj) [92] 

7/8 TeV 

o-84;°i2 

(0.87,0.13,0.0) 

hjj TT (2j) [92] 

7/8 TeV 

0-95;°;f8 

(.17, .83,0.0) 

Vh —>■ TT [92] 

7/8 TeV 

u.o/_o_88 

(0.0,1.0,0.0) 

h WW (0,lj) [92] 

7/8 TeV 

n 77+0.23 

U. 1 /_o .21 

(0.83,0.17,0.0) 

h WW (2j) [92] 

7/8 TeV 

0-62;°;f8 

(0.17,0.83,0.0) 

Vh WW [92] 

7/8 TeV 


(0.0,1.00,0.0) 

ZZ [92] 

7/8 TeV 

0.88;+27 

(0.9,0.1,0.0) 

h^ZZ (2j) [92] 

7/8 TeV 

1 cc+0.95 
-L.OO_o 66 

(0.58,0.42,0.0) 

(ggH) [94] 

7/8 TeV 

1 1 O+0.37 
-‘-•-*-^-0.32 

(1.0,0.0,0.0) 

/r ^ 77 (VBF) [94] 

7/8 TeV 

-L-00_o.68 

(0.0,1.0,0.0) 

^ 77 (VH) [94] 

7/8 TeV 

-o-ie+i® 

(0.0,1.0,0.0) 

ft. —77 (ttH) [94] 

7/8 TeV 

2.69;+8] 

(0.0,0.0,1.0) 


Table 6: CMS Higgs analyses used in constructing coupling fits. The best-fit signal strength 
modifier is denoted by fi with corresponding ±lcr errors. The relative contributions e are 
reported for production initiated by gluons via gluon fusion (GGH), weak gauge bosons 
via vector boson fusion or vector associated production (VBF/VH), and top quarks via ti 
associated production (TTH). 


We construct a combined likelihood from the product of all single-channel likelihoods, 

= (A.3) 

i 

This approach does not take into account correlations among systematic uncertainties in 
different Higgs searches, as such information is not publicly available. However, this is a 
reasonable approximation since uncertainties in Higgs measurements are not yet dominated 
by systematics. We are often interested in treating some inputs as nuisance parameters 9, in 
which case the combined likelihood may be expressed as a function of both /r and 9. 

We construct coupling fits using the profile likelihood approach [62]. In this approach, 
the best-fit signal strength modifier /I and corresponding uncertainty A/t of the combined 
likelihood are calculated using the likelihood ratio X{fi) = 9)/C{jl, 9). This is the ratio of 

a likelihood function with nuisance parameters 9 optimized for a given value of /r to a likelihood 
function where jl and 9 are optimized simultaneously. Optimizing nuisance parameters 9 for 
a given value of [i amounts to profiling these nuisance parameters. Given this likelihood ratio, 
the uncertainty A/i is computed using the test statistic —21nA(/r), which converges to a 
distribution in one degree of freedom as the data sample size is taken to be large. 
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B The gg ^ H/A ^ tt Differential Cross Section 


The kinematics of spin averaged two to two on-shell scattering processes 12 —)• 34 are com¬ 
pletely specified by the masses of the particles, energy-momentum conservation, and the 
Mandelstam variables s = {pi+P2)‘^ and t = (pi— Ps)^- In order to characterize the invariant 
phase space distribution for these processes it is useful to define a shifted dimensionless version 
of the Mandelstam variable t that has a flat metric with respect to the squared amplitude. 
For mi = m2 = 0 and m3 = m 4 = m this variable is 

w = 1-\—— m^'j = (3 cos 9 = tanh ^ ^ (B.l) 


where /3 = — 4m^/s is the velocity of either final state particle in the center of mass 

frame, cos 9 is the cosine of the center of mass scattering angle, and y = ia,ohr^{pz/E) are 
the individual rapidities of the final state particles in any longitudinal frame. The scattering 
variable lies in the range —/3 < w < (3. At fixed s the differential cross section with respect 
to w and t and cos 6 are related by 


da s da 1 da 

dw 2 dt {3 dcosO 


(B.2) 


Experimentally the variable w is more robust than cos0 near threshold since /? —)• 0 and 
cos 9 becomes ill-defined at threshold, but -oj —)■ 1 in this limit. In addition, s da/dw is 
proportional to the dimensionless scattering amplitude squared everywhere in phase space. 
So the dimensionless coordinates w and = s/(4m^) provide a flat metric for the prob¬ 
ability density in phase space, with the physical phase space region given by > I and 

— Y^I — 1/x^ < ro < s/l — I/x^. Central scattering corresponds to ro = 0. The phase space 
volume vanishes at threshold, = I, ro = 0. These coordinates for two to two scattering are 
the analogs of Dalitz coordinates for three-body decay. 

For heavy scalar and pseudo-scalar Higgs bosons with masses greater than twice the top 
quark mass, mu, aiA > 2mt, the scattering of two gluons to a top and anti-top quark, gg —>■ tt, 
receives contributions both from QCD interactions, as well as from s-channel gluon fusion 
production and decay of H and A. The differential cross section including all these processes 
is given by [39] 


da . 


Tvra^ 


2 r 


f{s/Amt,w)QCD + f{sl4m‘f, sl4ml,m\ilm\, ro)H-QCD 


-F f{s/4m^, s/4ml, m^mf, ro)A-QCD + j{s/4m\, sl4ml,m\lm\, w) 


H 


-F !{sl4m\,sl4ml,m\lm\,w)K 


(B.3) 
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where 


f{xt ,^)qcd = 
f{xt,xl,r^,w)u-QCD = 
/(a:?,a;b,r^,ro)A-QCD = 

f{xt,xl,r‘^,w)n = 

f{xt,xl,r‘^,ru)A = 


7 + 9w‘^ 
7(1 — w"^) 

24 
7(1 — 

6 

7(1 — wY 


1 2 

1 

2 ^ 

1 + w^ xf 

B 

1 

_I4e 

[4xf — r2 — 4iXfjH(xt)]N(Xf,x^) 

V xfJ 

(4xf — r2)2 + 16x^71^ (xt) 


Re 


72x2 


72x2 

112 


1 

1 - ^ 

xt 


[4x^ - r2 - 4ix^7A(xt)]P(x^, xf) 
(4x2 — r'^Y + 16xj7^(xt) 

[4x2 — Y — Aixl'^H{xt)]N{xl^xl) 


(4x2 — r2)2 + 16xj7|f(xt) 


l 2 


[4x2 _ ^2 _ 3-2^ 


(4x2 — r 2)2 + 16x^7^(xt) 


l 2 


(B.4) 


The first term for the QCD scattering amplitude squared function is unity for central scat¬ 
tering both at and well above threshold where /(1,0)qcd = /(c>o,0)qcd = 1- The second 
and third terms arise from interference between the heavy scalar and pseudo-scalar Higgs 
s-channel amplitudes and the spin and color singlet component of the QCD amplitude. Near 
threshold the scalar Higgs amplitude interferes with the P-wave component of the QCD am¬ 
plitude, while the pseudo-scalar amplitude interferes with the S'-wave component. These 
terms include one-loop functions for the gluon fusion production amplitudes generated by 
top and bottom quark loops 


N{xYxl) 

P{xt,xl) 


\/2 Gp 

47r^ 


f=t,b 



\/2 Gp 


A A 

mtKt y 

f=t,b 


In 


1 + ^1 _ 1/^2 
1 - ^1 - l/x 2 


1 + ^1 - l/x 2 
1-^1- l/xj 

. V 

— ^7^ 



2 


(B.5) 


The in terms come from absorptive branch cuts in the one-loop amplitudes. For the top quark 
loop contributions to the full amplitude, these absorptive pieces represent QCD production 
of an on-shell top and anti-top quark with final state re-scattering through intermediate 
s-channel heavy Higgs bosons. For niH = rriA mt the absorptive contributions to the 
interference of the heavy scalar and pseudo-scalar Higgs bosons with QCD are equal in both 
magnitude and sign for all center of mass scattering energies. Away from the heavy Higgs 
resonances, interference of the QCD and non-absorptive parts of the heavy Higgs amplitudes 
are suppressed compared with the QCD amplitude squared by 0{XtXf /4:n‘^) where f = t,b. 
The functional dependence of the interference terms on the angular scattering variable zu 
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Figure 21: Isocontours for the differential cross sections d?a/dvjdmii in QCD (red dotted 
lines) and for the interference terms in QCD+H (black solid lines) at 14 TeV with = 1. 
The left (right) panel shows the results for niH = 400 (600) GeV. 


is universal and independent of the center of mass scattering energy. The real part of the 
interference terms changes sign across the resonances, leading to a distinctive excess below, 
and deficit above, the resonances [39]. The fourth and fifth terms arise from the square of 
the amplitudes for s-channel production and decay of the heavy Higgs bosons. Away from 
the heavy Higgs resonances, the non-absorptive parts of the heavy Higgs amplitudes squared 
are suppressed compared with the QCD amplitude squared by 0(A|Aj/167r^) where f = t,b. 
The running dimensionless widths of the heavy Higgs bosons are 


7x(x) 


rx(^) 


(B.6) 


where Tx = r(A —)■ All) are the running total widths for X = H, A. These terms in the 
s-channel heavy Higgs propagators represent absorptive hnal state re-scattering of the heavy 
Higgs bosons through all intermediate on-shell states that contribute to the Higgs boson decay 
widths. 

The scalar and pseudo-scalar Higgs amplitudes do not interfere in the spin averaged top 
and anti-top quark production differential cross section. Interference could arise first in a 
spin averaged production differential cross section with at least two gluons radiated from the 
final state. It also arises in the spin averaged differential cross section of the full phase of the 
top and anti-top quark decay products. 

In Fig. 21 we plot contours of the differential gg —)• tt cross section d'^cj/dwdmu in QCD 
and of the interference contribution in QCD with a heavy scalar Higgs boson. We see that 
although the angular scattering variable offers some distinctive discrimination, there is no 
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region where magnitude of this discrimination is able to effectively overcome the prodigious 
QCD background. 


C Top Quark Reconstruction 


In the H/A ^ ti, hbH/A —)• bbtt, and charged Higgs analyses, we reconstruct the W bosons 
and top quarks with the following algorithm. 

The hadronically-decaying W boson is reconstructed using the non-6-tagged jets in the 
events. We choose the pair of the jets jij 2 which minimizes — mwhl: where rrij^^j^ is 

the invariant mass of the dijet system and m\Yh = 77.5 GeV.® The reconstructed hadronic 
decaying W 4-momentum is rescaled by a small correction factor of mw/mwhi where rrqy = 
80.4 GeV. 

To reconstruct the leptonically-decaying top quark, we solve for the 4-momentum of the 
neutrino in the final state using neutrino and W mass shell conditions. The solution for the 
z-component of the neutrino momentum is 




A{mw)P£L ± Ei^A [mwf - 


2p- 


iT 


(C.l) 


where A {my/) = When there are two real solutions, the mass-shell condition 

of the top quark may be used to select the “true” solution. Due to detector smearing effects 
and the finite width, there may be no real solution for For such events with negative 
discriminant, we enforce 

A {mw) = (C.2) 


and look for a modified I^t' which minimizes H is clear that when ^ ij' is 

parallel to I^Ti the discriminant must be positive. So we can expand 


= x'^ £T + (1 + 


(C.3) 


Then we have a constrained minimum value problem, which can be solved with Lagrange 
multipliers. We obtain a cubic surface and with a unique real solution 


y = - 


X = — 


(P - {5m^ -h 4c) d -h (m^ - 4c) ^ 


(P + (m^ - Ac) d+ {myy - 4c) ^ 

(P + (m^ -I- 8 c) d -I- (m^ - 4c) ^ 
12ad 


y 


(C.4) 

(C.5) 


°mwh is the center of the distribution obtained from the simulation of a pure sample of hadronic-VK 

events. 
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where 


a = 





(G.6) 

b = 





(G.7) 

c = 

'f IT ■ 




(G.8) 

d = 

|216a6?TT,^ -|- (jn^ -|- 8c) — ‘IQtyv^c - 

- 8c^) 



+ 

12\/3mvy 

WSabui^ -|- -\- 2c) 

- 16c) 

1/2 


X 

(a6-c2)V2}V3_ 



(G.9) 


We show the correction to the missing transverse momentum in figure 22. It is evident that 



S^/^T 

Figure 22: The correction to the missing transverse momentum used to obtain a real solution 
to the mass shell equation for the leptonically-decaying W boson. 


for most of the events with negative discriminant we only need to shift the missing transverse 
momentum by a factor of order 1%. 

To reconstruct the top quarks in the final state, we try all of the combinations {Whbh){Wibi), 
where Wh and Wi are the reconstructed hadronic and leptonic W bosons and bh and bi are 
6-tagged jets in the event. When there are two different real solutions of the neutrino longi¬ 
tudinal momentum, both of them are used. We select the combination which minimizes 


X = 


{mwhbh - , {rnwth - ^t) 


+ ■ 


crt 


crt 


(C.IO) 


where mt = 173.2 GeV is the pole mass of the top quark, ah = 50 GeV and ag = 25 GeV. 
To check the reconstruction efficiency, we compare the reconstructed top quark 4-momenta 
with the real parton-level momenta in the corresponding event. We calculate the ratio of 
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the modulus of the reconstructed top quark 3-momentum in the corresponding parton-level 
top quark rest frame to the energy of the parton-level top quark in the laboratory frame, 
result is shown in fignre 23 (note the logarithmic z-axis). It is clear that 
most of the reconstructed top quarks fall in the 6p/E < 0.15 region, meaning the top quarks 
in the events are well-reconstructed. 



Figure 23: The 2-dimensional distribution of the reconstruction efficiency of the hadronic 
and leptonic decaying top quarks in the events. 
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